Low-energy excitations of spin-Peierls chains with modified bond-impurities 



D. Augier", J. Riera'' and D. Poilblanc'^ 
Lahoratoire de Physique Quantique & Unite Mixte de Recherche 5626, Universite P. Sabatier, 31062 Toulouse, France 
^ Instituto de Fisica Rosario, Consejo Nacional de Investigaciones Cientificas y Tecnicas y Departamento de Fisica, 
Universidad Nacional de Rosario, Avenida Pellegrini 250, 2000-Rosario, Argentina 

(February 1, 2008) 



The introduction of modified bond-defects in spin-Peierls systems is investigated in a model of 
antiferromagnetic Heisenberg spin chains coupled to adiabatic phonons. Generically, new low-energy 
magnetic or non-magnetic excitations appear below the bulk spin gap energy. When two adjacent 
bonds are modified, these excitations can be interpreted in terms of bound states of a soliton with 
the localized spin-1/2 located on the impurity site. It is shown that the confining potential occurs 
even in the case of isolated chains. 

PACS numbers: 75.10 Jm, 75.40.Mg, 75.50.Ee, 64.70.Kb 



I. INTRODUCTION 

Impurity doping in quasi-one-dimensional spin-Peierls 
systems has recently sparkled renewed both experimen- 
tal and theoretical attention in this field. Experimen- 
tal studies in CuGeOs include so far in-chain magnetic 
(spin-1 Ni) or non-magnetic (Zn, Mg) impurity doping 
and ofF-chain (Si substituting Ge) doping. A rapid re- 
duction of the spin-Peierls temperature has been jgeneri- 
cally observed as impurity concentration increasead. Fur- 
thermore, impurity doping favors a new phase at low 
temperature, which has been characterized as a three- 
dimensional antifcrxomagnetic (AF) ordering by specific 
heat measuremftJitsEI, neutron scatteringB or nuclear mag- 
netic resonancea. More surprisingly, the AF phase has 
been shown to coexist with, ithc spin-Peierls phase for 
low impurity concentrationsHi3. For a theoretical under- 
standing of these effects, the starting point is the real- 
ization that the elementary excitations of these quasi- 
one-dimensional systems are objects which appear on 
chains and which have beeii^haracterized as topolog- 
ical defects called solitona3Tl3. These solitons do not 
interact with non-magnetic impurities in a strictly one- 
dimensional (ID) model, but are bound to spin-1 im- 
purities. However, impurities and solitons are confined 
into boimd states by the three-dimensional character of 
phononsEJ. 

Most theoretical studies have so far been concentrated 
on the effects of in-chain impurities in which a Cu ion 
is substituted by another ion. In this paper, we theo- 
retically study a new class of impurity systems, namely 
chains where some bonds have been modified. This could 
physically correspond to the substitution of a Cu ion by 
another spin-1/2 ion thus changing the values of the ex- 
change couplings connecting the impurity site with its 
two nearest neighbors (NN) sites. If next nearest neigh- 
bors (NNN) couplings are included in the model, they 
would be modified as well. The case of the spin-1/2 im- 
purity could also be thought as coming from defects in the 



crystal structure (structural disorder, radiation damage, 
grain boundary in polycrystalline samples, etc.) which 
change locally the values of the couplings. These modi- 
fied couplings are then described by an effective spin-1/2 
impurity (see e.g. Ref. ^2|). Similar kind of defects could 
also be due to the substitution of the out-of-chain Ge ions 
by Si (see e.g., Ref. H) then changing only one NN and 
two NNN bonds ("bond-centered impurity"). Since, as 
it was pointed out above, the elementary excitations in 
these quasi- ID systems are objects which exist already 
in strictly ID systems, the first step is to characterize 
these excitations by using purely ID models and this is 
the purpose of the present work. A second step would 
be then to understand how these excitations and their 
interactions are modified by interchain (magnetic and 
elastic) interactions present in CuGeOa. The first study 
of the effects of bond impurities has been done on iso- 
lated Heisenberg chains in Ref. |l^ using essentially field- 
theoretical techniques. We will compare our results with 
the ones obtained in that study in Section ^ For spin- 
Peierls systems, analytical studies on the same model we 



will consider (see Eq. I 
bosonization technique: 



L- below) have been done using 
llil. We think, however that in 
those studies there is not a proper understanding of the 
elementary interactions at a microscopic level. 

In practice, we introduce these bond impurities in our 
model by redefining the exchange coupling constants on 
two bonds connected at the impurity site (spin-1/2 or 
"site-centered impurity" ) or a single bond (bond-centered 
impurity) as Jimp = x J. In the first case, from a purely 
theoretical point of view, this allows us to interpolate 
continuously between a periodic chain with even number 
of sites and an open chain with odd number of sites plus 
a spin disconnected to the chain. We find that there are 
also some interesting features in the region a; > 1. By 
assuming that lattice distortions are adiabatic, we are 
thus lead to the following Hamiltonian : 
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where i (j, I) indicate bulk bonds (impurity bonds, all 
bonds, respectively). We assume for simplicity that 
the spring constants K are not modified by the defects. 
Moreover, we assume also that the two NNN couplings 
at the impurity site also involve the same frustration pa- 
rameter a, i.e. J/jjjp = a Jimp- The variable x will be 
extended to negative values in which case the impurity 
couplings are ferromagnetic, thus including the possibil- 
ity of studying a spin-3/2 impurity when two adjacent 
bonds are modified (e.g. a Co ion replacing a Cu ion in 
the chain, see Ref. pq) or a spin-1 Ni impurity if only one 
bond is changed (see Ref. 16). We have solved this Hamil- 
tonian numerically on finite chains by the self-consistent 
procedure described in Ref. |l^ using Exact Diagonaliza- 
tion (ED) and Quantum Monte Carlo (QMC). The rele- 
vance of this kind of calculations to describe experimental 
results on CuGeOa, in particular the use of the adiabatic 
approximation which somehow interpolates between the 
full quantum spin-lattice coupling and the simplest mod- 
els with fixed dimerig atipn , was emphasized in a number 
of previous studies£3'oEj la. particular, recent calcula- 
tions within this formalismEj using parameters realistic 
for CuGeOa, have given results that describe Raman ex- 
periments not only at a qualitative but at a quantitative 
level, in spite of the fact that for these experimepis dy- 
namical lattice effects are in principle important. Eil 



II. CHAINS WITH TWO ADJACENT MODIFIED 
BONDS 

We start our study with the case of a "site-centered" 
impurity. A frustrating exchange Jj'j^^p = aJimp = ctxJ 
is also included between the impurity site and its two 
next nearest neighbors. We first consider even chains 
in the region < x < 1. In the case of a = 0, the 
X — point is equivalent to an open odd chain with 
a disconnected spin- 1/2. It has been previously shown 
that, within this ID model, no binding occurs betweea. 
the spin- 1/2 soliton released in the chain and its edgestll 
(which here correspond to the location of the impurities). 
On the other hand, for x = 1, we recover the pure system 
and there are of course no magnetoelastic solitons in the 
ground state (GS). One can then consider that the soliton 
is tightly bound to the spin-1/2 impurity, forming a spin 
singlet. As x is reduced from 1 to 0, this bound state 
is weakened and the soliton eventually moves away from 
the impurity site. This behavior can be qualitatively seen 
from the QMC data at T 0.05J for L ^ 40, K ^ 
1, as shown in Fig. 

0- In (a), X — 0.9, the distortion 
pattern is almost indistinguishable from that of a pure 
chain with periodic boundary conditions. In (b), x = 0.3, 



the soliton has only moved one lattice spacing away from 
the impurity. In contrast, as shown in (c) corresponding 
to X — 0.1, the bound state has disappeared and the 
soliton moves freely in a region around the center of the 
chain segment. 
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FIG. 1. QMC results at T = 0.05 J for the distortion 5, on 
a 40 site chain with PBC, K = 2, and for various values of 
x; (a) x=0.9, (b) x=0.3, (c) x=0.1, (d) x=-2.5. The spin-1/2 
impurity is located at site number 40 at the extreme right of 
the plots. 

To put these statements on a more quantitative basis, 
let us start by analyzing the energy per site as a function 
of X, e{x) = E{x)/L. Since for x = there is no soliton- 
impurity binding, a value of the energy lower than e{x = 
0) could be considered an indication of the presence of a 
bound state in the system. The QMC results for L = 40, 
K = 1,2 and 3, obtained for T — 0.05 which can be taken 
as the ground state energies for the system considered 
here, except very close to a; = as discussed below, are 
shown in Fig. ||. For all the values of K studied here it 
can be seen that the e{x) is maximum at a; = 0. This 
might indicate that there is an impurity-soliton bound 
state for cc > and that the binding energy decreases 
as x is reduced from x — 1, consistently with Fig. ^a), 
(b). The fact that e{x) is lower than e(0) for x < 
can be interpreted by assuming that an effective spin-3/2 
impurity starts to form and hence the soliton is again 
attracted to this impurity. For the particular value of 
K — 2,we have studied the size dependence of the results 
by computing the energy for L = 20 and L = 80 in 
addition to L — 40. The energy is always maximum at 
X = and, as expected, it tends to become constant as 
L increases. The insert shows the energy of the lowest 
state in the = 1 sector together with the energy of 
the excited state containing one extra soliton-antisoliton 
(s — s) pair. It can be seen that for a; < the ground 
state of the system has 5* = 1. We will discuss these 
features later. 
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FIG. 2. QMC results for the energy per site as a function of 
a; on a 40 sites chain (except otherwise stated) with PBC and 
one spinl/2 impurity, T = 0.05 and various values of K as 
indicated in the figure. The insert shows the first two = 1 
states for K = 2 and L = 40. 

The deconfincmcnt of the impurity-soliton bound state 
when a; ^ (a; > 0) is a rather subtle issue. The above 
results obtained at finite temperature T = 0.05J and 
X — 0.1 shown in Fig. [^(c) indicate that there is no 
soliton-impurity binding in that case. However, for a 
temperature larger than the (zero temperature) binding 
energy, one expects that, due to thermal fluctuations, the 
soliton can escape from the small confining potential. In 
order to confirm this scenario and to determine if the 
above bound state survive in the thermodynamic limit at 
zero temperature, we have computed the soliton-impurity 
(T = 0) binding energy which is rigorously defined as 

Eb{L) = {Ei„jp{L) — Eo{L)) — Cimp ~ Csoi, 

where e^oi = limL^oo[£^o(i + 1) - E^{L + 1)], EoiL) 
{Einip{L)) is the GS energy of a pure (bond-modified) L- 
site chain (L = 2p even), E^{L + 1) = (EoiL) + Eq{L + 
2))/2, and ei„,p = limL^oo[i^imp(^ + 1) - + 1)]. 

These quantities have been computed within an ED 
treatmendiil. An impurity-soliton bound state corre- 
sponds to Eb{oo) < 0. 

Some results for Eb{L) are shown in Fig. ^. For x ~ 0, 
by using extrapolations of the form a + b/L exp(— L/Lq), 
we have obtained a vanishing binding energy (within an 
estimated error of ~ O.OOIJ) as it should be, as soli- 
tons and non-magnetic impurities do not—bind in spin- 
Peierls chains including lattice relaxationlJ. For x = 1, 
the binding energy is the energy necessary to free one spin 
which is part of a singlet, i.e. the spin gap. We obtained 
Eb = — 0.39J = — Agpin, "vebich impHes that there is no 
s — s binding in this modc]E3. For x = 0.1, the situation, 
as depicted also in Fig. || is more intricate. In this case, 
there is a cross-over between the free and bound regimes 



near L = 14 as signaled by a change of the curvature of 
the scaling form of the binding energy. The exponential 
fit is of course no more valid but one can crudely esti- 
mate that in the bulk limit Eb < — O.OIJ. This result 
suggests that, for fixed positive x, binding sets up only 
for sufficiently large chain lengths. Alternatively, this im- 
plies the existence of a finite critical value of a;, XcriL), 
below which there is no impurity-soliton bound state for 
a given chain size L. This critical value Xcr{L) has been 
evaluated by ED for various lattice sizes and results are 
shown in Fig. ^. Note that this crossover disappears for 
X > 0.3 i.e. the soliton is bound for all sizes we consid- 
ered and the extrapolated binding energy is clearly non 
zero (negative). The extrapolation of the Xcr{L) curve 
when L — > oo is consistent with a finite binding energy 
in the thermodynamic limit for any finite value of x. 
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FIG. 3. Impurity-soliton bound state energy as a function 
of the inverse lattice size 1/L for a = K = J and x=0.0 
(o) or x=0.1 (□) obtained by ED. An exponential fit is also 
plotted for the case a; = 0. The arrow points out the crossover 
between two scaling behaviors. 

At finite temperature, a similar cross-over should also 
appear. However, since for sufficiently large chains the 
thermal length becomes the only relevant length scale, 
one expects that XcriL) extrapolates, when L — > oo, to 
a (small) finite value consistently with the above QMC 
results of Fig. |l](c) and in contrast with the T=0 case. 
The schematic behavior of Xcr{L) at finite temperature 
is shown in Fig. ^ for comparison. 

Let us now consider negative values of x, i.e. the case 
of a ferromagnetic bond. The limit |a::| >> 1 is simple to 
understand. In this case, a spin-3/2 effective spin forms 
around the impurity leaving a spin- 1/2 in the rest of the 
even chain. This behavior is qualitatively shown in (d) for 
X = —2.5. Since, as we shall discuss more quantitatively 
later on, the soliton is weakly antiferromagnetically cou- 
pled to the spin-3/2 impurity, we expect the GS to be 
degenerate. On the contrary, for |a;| << 1, the small 
ferromagnetic couplings connected to the impurity site 
produce an effective interaction that weakly binds the 
localized spin-1/2 to the spin-1/2 soliton into a triplet 
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bound state. 
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FIG. 4. Evaluation of critical parameter Xcr as a function 
of the inverse lattice size 1/L by ED. Parameters such as 
a = and K = J have been used. An hypothetic curve at 
finite T is also plotted. 

In Fig. I we show the zero temperature binding energy 
in the thermodynamic limit as a function of x together 
with the estimate of the error (resulting from the finite 
size scaling procedure). These results are consistent with 
the above qualitative discussion and confirm the exis- 
tence of a soliton-impurity bound state in the bulk limit, 
except for x = 0. In the x < part, a spin-3/2 is located 
on the impurity site and its two neighbors, and a finite 
binding is also observed. In the x > 1 region, the soliton 
is very strongly bound to the impurity, and it becomes 
impossible to deconfine the two objects. In other words, 
trying to pull out the two components of the singlet pair 
from each other will result into the spontaneous creation 
of a soliton-antisoliton pair out of the "vacuum" . Con- 
sequently, the binding energy, as it is defined, becomes 
identical to the spin gap i.e. the energy to create such 
s — s pair. 




To complete the overall picture, we have computed 
the singlet-triplet spin gap defined as = £'(5' — 
1, {(5}i) - E{S^ = 0, {<5}o) where {S}o, {S}i are the set of 
bond distortion determined for S — and 5 = 1 respec- 
tively. It is striking that the overall behavior of vs x 
is totally similar to that of the binding energy shown in 
Fig. ||. This strongly suggests that around a; = the low- 
est singlet {x < 0) or triplet {x > 0) excitation is made 
of a deconfined soliton. When a; > 1, the spin gap is con- 
stant as the soliton remains bound to the impurity, and 
bulk s — s excitations are energetically more favorable. 
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FIG. 6. Energy difference between the S = 1 and S = 
sectors obtained by ED in the bulk limit as a function of x 
for a — and K — J. 

It is instructive to relateJiie above results to the anal- 
ysis of Eggert and AfHeckllj for a spin-1/2 Heisenberg 
chain (without coupling to the lattice) in the presence of 
impurities. In this work, they show that when a; > 0, in 
the renormalization group language, the open chain will 
be unstable and ultimately flow to the stable periodic 
chain with the impurity site included ("healing" effect). 
This is consistent with our result that there is a bound 
state for all a; > in the bulk limit. Now, for a; < 0, the 
system would have a marginal flow (as i — > oo) towards 
the open chain with a decoupled spin, i.e. no impurity 
soliton-binding in the bulk limit. This result is the oppo- 
site to our result shown in Fig. ||. This indicates that the 
Luttinger liquid approach of Ref. |l3| cannot be directly 
extended to the case when a coupling to the lattice is 
present. 



III. CHAINS WITH A SINGLE MODIFIED BOND 



FIG. 5. Impurity-soliton binding energy obtained by ED 
in the bulk limit as a function of x. Parameters a — and 
K = J have been used. 



We now turn to the case where a single bond is mod- 
ified, i.e., on one NN bond, J Jimp: a-nd its two over- 
lapping NNN bonds, aJ Jimp- This configuration 
corresponds to a "bond-centered" rather than a "site- 
centered" impurity. 

In the present situation, even in the case a; = 0, the 
number of sites in the chain remains even and so there is 
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no extra spin carrying the soliton as in the previous case. 
Hence, QMC results for the distortion Si, ior < x < 1, 
for K = 2 and a = look very similar to those depicted 
m Fig.|] (a). As it was mentioned in the Introduction, 
in the limit x ^ ~oo, the impurity bond leads to an 
effective 5 = 1 impurity. On the other hand, for x oo, 
a strong dimer is formed at the modified bond and it will 
be relatively decoupled from the remainder of the chain. 
However, for intermediate values of x, it is possible that 
a s — s pair can be bound to the impurity bond and, in 
fact, QMC results for Si indicate that this may be the 
case for a; < 0. 
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FIG. 7. QMC results for the energy as a function of x on 
a 40 sites (except otherwise stated) chain with PBC and one 
modified bond, T = 0.05 and various values of K as indicated 
in the figure. The insert makes a comparison of the energies 
of the lowest states in the 5*^ = and — 1 sectors for 
L = 40, and K ^ 2. 

If the energies of large systems are examined (FigJ^), 
it can be seen that the highest energy is again located at 
X = 0. Notice that for x < the system is frustrated and 
classically a higher energy than for a; = is expected. 
For all X, the GS is a singlet and the 5=1 states are 
well separated from it, as the insert of Fig. ^ shows. Also 
as this insert indicates, there is a crossover at a; = 1 be- 
tween two — (in fact 5 = 0) states. The ground 
state for a; > 1 corresponds to the above mentioned rela- 
tively decoupled dimer which leads to an energy decreas- 
ing roughly linearly with a:. It is then clear that the local 
"impurity" potential lifts the degeneracy of the two-fold 
degenerate GS occurring exactly at a: = 1 immediately as 
one moves away from this particular point. Indeed, the 
two states crossing at a: = 1 exhibit opposite bulk dimer- 
ization far away from the impurity site. This can be 
qualitatively understood from the fact that, for a; < 0, 
a singlet is formed around the impurity which involves 
three bonds weakly connected to the rest (a triplet on 
the center bond bound to two solitons) compared to only 



one when a; > 1. In the limit a; — > — oo one recovers 
the case of a spin-1 impurity (corresponding to a Ni — > 
Cu substitution) previously studied.E2l Precisely in Ref. 
^ it was emphasized the role of the three-site subsystem 
composed by the S = 1 site and its two S = 1/2 NN sites. 

We have also computed the related spin gap for this 
type of impurity. Results obtained by both ED and QMC 
are shown in Fig. ^j. The extrapolation to the bulk limit 
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FIG. 8. ED and QMC results for the spin triplet gap as 
a function of x on periodic chains and one modified bond for 
K = 2. The extrapolation to the bulk is also shown in the 
interval —2 < x < 2. 

using the law a exp(— L/Lq)/-^'', is also shown. It can 
be seen that the gap is roughly constant (within error 
bars) in the interval —2 < a: < 2, in striking contrast with 
the behavior of the site-centered impurity as discussed 
in the previous Section. It is also interesting to notice 
that two different triplet excitations also cross exactly 
at a; = 1. Since the energy separation from the GS is, 
to a good accuracy, independent of x and equal to the 
bulk spin gap, these states can be interpreted as a bulk 
(deconfincd) s — s pair excitation. 

IV. SUMMARY AND CONCLUSIONS 

Based on the previous QMC and ED results, a 
schematic representation of the behavior of the low- 
energy states with x can be drawn. A summary of the 
structure of the low-energy spectrum in the case of a 
"site-centered" impurity is schematically represented in 
Fig. ||. For X > 1, the GS includes an impurity region 
consisting of two adjacent strongly dimerized bonds car- 
rying a spin-1/2 and separating two dimerized regions 
with opposite q = n lattice order parameters. On an even 
chain, this defect binds a soliton. Since, on a closed chain, 
the soliton can be on either side of the defect the GS is 
then two-fold degenerate as shown in Fig. a The lowest 
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excitation corresponds to a deconfined soliton-antisoliton 
pair decoupled from the impurity (which could be in sin- 
glet or triplet degenerate states) as represented by the 
dotted line in Fig. ||. However, for < x < 1, the low- 
est excitation is a quite different state involving a spin 
flip of the soliton spin linked to the spin-1/2 defect at 
the impurity site. This qualitative change manifests it- 
self as a level crossing of the first excitations at x = 1. 
At exactly a; = 0, we expect additional level crossings 
of both GS and excited states as shown in Fig. |[ In- 
deed, in this case, the first excited state in a chain with 
an even number of sites (including the impurity site) is 
highly degenerate and corresponds to the 3 possible spin 
configurations (5 — 0,1,2) of the three free spin-1/2 of 
the two mobile solitons and the localized spin-1/2. Arbi- 
trary small ferromagnetic or antiferromagnetic couplings 
on the two impurity bonds are expected to lift this de- 
generacy into several levels as seen in Fig. ^in agreement 
with the previous numerical calculations. 




FIG. 9. Schematic representation of the low-energy spec- 
trum of an even chain with a "site-centered" defect as a func- 
tion of X. The spin sectors as well as the nature of the states 
are indicated in the figure. 

A similar plot of the low-energy spectrum is shown in 
Fig. ^ in the case of a "bond-centered" impurity. Con- 
trary to the previous case, such a defect favors one of the 
two bulk dimerization patterns so that the GS is non- 
degenerate for all X, except at x = 1. For x > 1 the 
dimerization pattern with a singlet located on the modi- 
fied (strongly dimcrized) bond is selected while, for x < 1, 
strong bonds occur next to the modified one. These two 
orthogonal GS naturally cross at a: = 1 as seen in Fig. |l^. 
For all x, the lowest triplet excitation corresponds to the 
creation of a bulk (deconfined) s — s pair. 

We finish this work by discussing the implications of 
the above results to experimental systems of doped spin- 
Peierls chains. Impurity doping should lead to important 
effects in magnetic properties (i.e. in inelastic neutron 
scattering, Raman scattering, etc..) when the spin ex- 
citation gap is reduced, < Agpin (0 < a; < 1) or 
when the GS carries a finite spin {x < 0). For defects 
centered on a site affecting the values of the exchange 



couplings on each side, low-energy magnetic excitations 
are predicted by our model as seen in Fig. ^. Such states 
seen in inelastic neutron scattering of Zn-doped CuGeOa 
materials (x — 0) should also survive for more general 
kinds of defects such as substitution of Cu by a different 
spin- 1/2 impurity. 
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FIG. 10. Schematic representation of the low-energy spec- 
trum of an even chain with a "bond-centered" impurity as a 
function of x. The spin sectors as well as the nature of the 
states are indicated in the figure. 

The substitution of a Ge ion by a Si one would pro- 
duce a decrease of the neighbor Cu-O-Cu angle and at the 
same time presumably a shortening of the corresponding 
Cu-Cu distance.El The outcome of these competing ef- 
fects on the value of the effective NN Cu-Cu exchange 
coupling is difficult to predict but it is quite likely that 
X = 1-1- Aa; with \Ax\ < 1. For both signs of A^; there is a 
common behavior: the impurity bond forces the chain to 
pick one of the two possible dimerization patterns. This 
selected pattern may or may not correspond to the one 
determined for the whole plane due to interchain mag- 
netic and elastic interactions. In the former case, there 
is a cost in energy which is proportional to the length 
of the "wrong" segment of the chain (assumed finite for 
a finite defect density). Then, this energy cost may be 
high enough to allow the formation of a soliton in an odd 
chain or a s — s pair on an even chain. The solitons will be 
bound to the defect and at this point the analysis follows 
the one pre*dously developed for the case of non-magnetic 
impurities. t2l Thus, we have for Si-doped CuGeOs the 
main behavior seen for in-chain non-magnetic impurities, 
nanaelK a transition for the spin-Peierls phase to an AF 
one.cTEl Moreover, in the vicinity of x = 1, we predict 
a new low-energy singlet excitation below the spin gap 
as seen in Fig. |l^. Such an excitation should not be 
affected by interchain interactions. It is tempting to re- 
late this low-energy singlet excitation to the ones seen in 
Raman scattering experiments-B However, a more quan- 
titative study is necessary to uniquely characterize this 
excitation. Besides, other changes detected in these ex- 
periments, like the decrease of the continuum intensity 
with increasing doping were rather attributed to an en- 
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larged NNN interaction along the chains, an effect which 
was not included in our study. 

ACKNOWLEDGMENTS 

We thank IDRIS (Orsay) for allocation of CPU time 
on the CRAY supercomputers and rs6000 workstations. 
J. R. also thanks SCRI and ACNS (Florida State Uni- 
versity) for using their computing facilities. 



D. Augier, D. Poilblanc, E. S0rensen and I. Affleck, Phys. 
Rev. B 58, 9110 (1998). 

P. Hansen, J. Riera, A. Delia, and E. Dagotto, Phys. Rev. 
B 58, 6258 (1999). 



^ M. Hase, I. Terasaki, Y. Sasago, K. Uchinokura and 
H. Obara, Phys. Rev. Lett. 71, 4059 (1993); K. M. Ko- 
jima, et al, Phys. Rev. Lett. 79, 503 (1997). 

2 B. Oseroff, S.-W. Cheong, B. Aktas, M. F. Hundley Z. Fisk 
and L. W. Rupp, Phys. Rev. Lett. 74, 1450 (1995). 

^ J. G. Lussier, S. M. Coad, D. F. McMorrow and D. McK. 
Paul, J. Phys. Cond. Matt. 7, L325 (1995). 

* J. P. Renard, K. Le Dang, P. Veillet, G. Dhalenne, 
A. Revcolevschi and L. P. Regnault, Europhys. Lett. 30, 
475 (1995). 

^ L. P. Regnault, J. P. Renard, G. Dhalenne and A. 
Revcolevschi, Europhys. Lett. 32, 579 (1995). 

''p. Lemmens, et at, Phys. Rev. B 55, 15076 (1997). 
M. Weiden, R. Hauptmann, W. Richter, C. Geibel, P. Hell- 
mann, M. Koppen, F. Steglich, M. Fischer, P. Lemmens, 

G. Giintherodt, A. Krimmel, and G. Nieva, Phys. Rev. B 
55, 15067 (1997). 

* B. S. Shastry and B. Sutherland, Phys. Rev. Lett. 47, 964 
(1981). 

^ T. Nakano and H. Fukuyama, J. Phys. Soc. Japan 49, 1679 
(1980). 

D. Khomskii, W. Geertsma and M. Mostovoy, Czech. J. 
Phys. 46, Suppl. S6, 3239 (1996). 

P. Hansen, D. Augier, J. Riera, and D. Poilblanc, Phys. 

Rev. B 59, 13 557 (1999). 

A. Smirnov et al, |cond-mat/9806347 . 
" S. Eggert and L Affleck, Phys. Rev. B 46, 10 866 (1992). 

H. Fukuyama, T. Tanimoto, and M. Saito, J. Phys. Soc. 
Japan 65, 1182 (1996). 

P. E. Anderson, J. Z. Liu, and R. N. Shehon, Phys. Rev. 
B 56, 11 014 (1998). 

S. Coad, J.-G. Lussier, D. F. McMorrow, and D. McK. 



Paul, J. Phys. 8, 6251 (1996); V. N. Glazkov, et al, cond- 



mat/9806325. 



""^A. E. Feiguin, J. Riera, A. Dobry and H. A. Ceccatto, Phys. 
Rev. B 56, 14 607 (1997). 

M. Fabrizio and R. Melin, Phys. Rev. B 56, 10 866 (1992). 
A. Dobry, P. Hansen, J. Riera, D. Augier, and D. Poilblanc, 
Phys. Rev. B 60, 4095 (1999). 

D. Augier, E. S0rensen, J. Riera, and D. Poilblanc, Phys. 
Rev. B 60, 1075 (1999). 

G. Els et al, Europhys. Lett. 43, 463 (1998). 



7 



